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(57) Abstract 

An active control method and apparatus for controlling vibration 
or sound wherein the computational bimicn to adaptivcly update the 
control filter within the update model is reduced. The apparatus includes 
means (not represented) for producing an input signal (36) for input 
to the update model (22) and the reduction block (26). means (not 
represented) for generating an error signal (41). and an output device 
(25) for canceling sound and/or vibration at a selected location. The 
method uses a gradient descent algorithm such as the filtered-x LMS 
algorithm and replaces the long filter model (23. 24) with a shortened 
length filter model (28, 30, 30'). Error signal infomnation together with 
the output infonnation from the shortened length filter model are used to 
update the coefficients of the control filter (32) according to an update 
calculation method. The shortened filter model only contains frequency 
response information at the L frequencies of interest. Because of the 
reduced amount of information contained in the shortened filter model, 
more updates per unit of time can be performed on the control filter(s) 
or a smaller or fewer processors or CPU's may be utilized. The method 
has applicability to feed -forward and feed-back type systems and also 
to systems using a waveform generator for supplying the input signal. 
The reduction method may be implemented on-line or off-line. 
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REDUCTION OF COMPUTATIONAL BURDEN OF ADAFTIVELY 
UPDATING CONTROL FILTER{S) IN ACTIVE SYSTEMS 

Field of the Invention 

5 This invention relates to the area of active control methods and 

apparatus for control of noise and vibration. Specifically, the invention 
relates to means and apparatus for adaptively updating one or more control 
filters used in active noise and/or vibration control systems. 

Background of the Invention 

IQ Active vibration isolation and sound isolation systems inject active 

forces and/or acoustic noise to effectively cancel noise and/or vibration at 
one or more defined locations. Hereafter, in the description, the plural will 
be used as meaning one or more. The noise and/or vibration emanates 
from a source of disturbance such as an engine. These active systems 

15 utilize microprocessors to calculate the appropriate response to eliminate, 
or at least to substantially minimize, mechanical vibration or ambient noise 
levels. Typical examples of active systems include systems which reduce 
noise and/or vibration experienced within an aircraft cabin or within an 
automobile passenger compartment. Generally, these active control 

20 systems are responsive to at least one input or reference signal and/or an 
error sensor signal. Typical devices for supplying these reference and 
error signals include microphones, accelerometers, force sensors, and the 
like. 

Examples of active systems are taught in US Pat. Nos. 4,677,676 and 
25 4,677,677 to Eriksson, the disclosures of which are hereby incorporated by 
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reference herein. Further examples include US Pat. Nos. 4,417,098, 
4,153,815 and 4,122,303 to ChapUn et al., 4,232,381 to Rennick et al., 4,083,433 
to Geohegan, Jr. et al., 4,878,188 to Ziegler, Jr., 4,562,589 and 4,473,906 to 
Waraaka et al., 5,170,433 to EUiott et al., 4,689,821 to Salikuddin et al., 
5,174,552 to Hodgson, and 5,133,527 to Chen et al. In these types of active 
systems, the algorithm used for control can be least mean square (LMS), 
Filtered-X LMS, Filtered-U, Newton's Raphson method, or other like 
gradient descent methods. For background, the Filtered-x LMS algorithm 
is described herein. 



10 



Filtered-x LMS Algorithm Background 



In active control systems, control filters such as finite impulse 
response (FIR) filters or infinite impulse response (IIR) filters are used. 
For example, a filtered-x system will include at least one control filter and 
at least one x-filter. The x-filter models the frequency response of the plant 
15 whereas the control filter transforms the input or reference signal into an 
output signal for driving an output device. In a single-input single-output 
(SISO) system illustrated in Fig. 1 utilizing the Filtered-x Least Means 
Square (LMS) algorithm, there is a feed-forward control filter AN,k and an 
x-filter Cm- The x-filter information and error signal ek information 
20 from an error sensor are used in an update calculation method utilizing an 
equation such as: Ak^i = Ak - ^ et Rk- The result of the update calculation 
method is used to adaptively update the coefficients or taps of the feed- 
forward control filter AN,k as indicated by the arrow. The Filtered-x LMS 
algorithm is fiilly described in "Adaptive Signal Processing" by B. Widrow 
25 and S. D. Steam, 1985. 

"Prior Art" On-Line Identification System 
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Fig. 2 illustrates an active c .rol system of the prior art which 
includes on-line identification of the x-filter model. This method is fully 
described in relation to Fig. 20 of US Patent No. 4,677,676 to Eriksson. 
Utilizing the numerals used in the *676 patent for reference, the Eriksson 
5 invention provides an active attenuation system and method for attenuating 
an undesirable acoustic wave by producing a canceling output acoustic 
wave from an output transducer 14. The key is that the invention models 
the feedback path 20, error path 56 and the output transducer 14 in an on- 
line fashion. In essence, the pre-training step previously required was 
10 eliminated, such that the active system may update itself for changes in the 
dynamic system due to temperature, etc. 

The acoustic system is modeled by adaptive filter model 40 provided 
by filters 12 and 22 and having a model input at 42. The model 40 outputs a 
correction signal at 46 to the output transducer 14 to introduce a correction 

15 signal to drive the error signal 44 to zero. The auxiliary noise source 140 
introduces noise into the adaptive filter model 40. The noise source is a low 
level random uncorrelated noise source. Model 142 models the error path 
56 on-line. A copy of S'E' is copied to 144 and 146 to compensate for the 
output transducer 14 and the error path 56. The error signal 44 is used 

20 with the output of the copy 144 and 146 at junctions 72 and 76 to provide 
update signals to the A and B filters. The novel featiire of this invention is 
that it requires no calibration, no pre-training of weights, and no start-up 
procedure as all updates are achieved on-line. 

"Prior ArT MIMO Systems 

25 In multiple-input multiple-output (MIMO) systems such as 

described in US Pat. No. 5,216,722 to Popovich and US Pat. No. 5,216,721 to 
Melton, the disclosures of which are hereby incorporated by reference 

3 
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herein, a feed-forward control filter is associated with each reference 
sensor/output device pair. Also, there is one x-filter associated with each 
output device/sensor pair. The coefficients or tap weights of each feed- 
forward control filter must be calculated and updated in real time. 

5 A typical MIMO system would have Na output devices, Ng error 

sensors, and Nr reference sensors. As a result, there will be (Na x Nr) 
control filters and (Na x Ns) x-Filters. For example, systems utilizing 16 
error sensors, 16 reference sensors, and 16 output devices, such as might be 
found in a typicaJ MIMO system, there will be 256 control filters and 256 x- 
10 filters. Each of the x-filter models may be developed on-line or off-line by 
pre-training and represents an estimate of the impulse response between 
each corresponding error sensor and output device. A matrix of 
information from the reference sensors is input into the x-filters. The x- 
filter output information is then used, along with the error sensor matrix 

15 information, to calculate the updated coefficients for each of the feed- 
forward control filters according to an update method. Agadn, this is 
known as adaptively updating the control filter tap weights. 

Typically, the x-filters are such that they model the plant over the 
entire frequency raxige of interest. Because of this, a large stream of 

20 frequency response information is developed to model the plant over this 
wide frequency range. These data streams may include hundreds of 
samples. As a result, the amount of information provided to the update 
calculation from the x-filters can be enormous. This results in the need for 
very fast, expensive microprocessors, or multiple microprocessors or 

25 CPU's to provide the required computing power to be able to update the 
control filter coefficients adequately. Because the computational burden for 
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MIMO systems is so large, typically the control filters are only updated in a 
round-robin fashion. 

In other words, not all the control filters are updated each update 
cycle in MIMO systems. For example, only two control filters may be 

5 updated each update cycle, and a different two control filters are updated 
the next cycle. As should be apparent, the computational burden becomes 
very large for MIMO systems, and as a result, adaptively updating of the 
feed-forward control filters cannot be accomplished very frequently. This 
causes both the convergence time and the level of noise/vibration 

10 cancellation achievable to suffer. 
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Summary of the Invention 

The present invention provides the method and apparatus for 
reducing the computational burden required to update control filters in an 
active control system by reducing the size of the filter models that are used 
5 for providing information to the update calculation method. For example, 
in x-filter systems, reducing the length, i.e., the size, of the x-filters allows 
a system to adapt to changes in the error path more effectively and quickly. 
In another aspect of the present invention, the reduction technique for 
reducing the size of the models takes place in an on-line fashion within a 
10 reduction block. This method and apparatus has particular applicability 
for systems having four or less discreet and dominant tones, i.e., tonal 
systems. However, the apparatus and method may also be used in systems 
having many tones. 



15 and therefore reduces the computational burden required of the 
microprocessor to perform the update calculation to update the control 
filters. These reduced or shortened x-filter models can have up to 25 times 
fewer values than prior art long x-filter models. This is especially true in 
tonal systems having a few discreet and dominant tones to be controlled. 

20 Fewer values means reduced number of calculations are required to 
adaptively update the control filters. Because of the shortened length, more 
updates per \mit of time can be achieved, thus providing better control of 
noise and/or vibration, or in the alternative, allowing the use of smaller and 
less expensive micro-processors or a smaller number of microprocessors. 

25 Of course all of the above making reference to filtered-x feed-forward 
systems have applicability to feed-back type systems also, as will be 



For x-filter systems, this method reduces the length of the x-filters 
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described later. Further, the present invention reduction method has 
appUcability in MIMO systems as well. 

The key to the present invention is that the long filter models are 
transformed into shortened or reduced filter models within the reduction 
5 block. The reduction method has particular applicability in tonal systems 
that contain four or less discreet and dominant tones. As a result of the 
smaller dimensional size, the calculations required to update the 
coefficients of the control filters are reduced correspondingly. For example, 
in x-filter systems, the long x-filter models are transformed into ^hort 

10 x-filter models which contain the frequency response information only at 
L firequencies of interest. In most aircraft control systems, there may only 
beL = 2toL = 4 dominant tones of interest. By extracting the frequency 
response estimates only at these L selected tones, and ignoring the other 
frequency response information, the length (size) of the x-filters can be 

15 dramatically reduced. Typical prior art frequency response models, i.e., 
long filter models may contain approximately 60 to 100 weights, whereas 
the reduced filter models may only require approximately 4 to 8 weights. 

The shortened filter models can be obtained by a relatively simple 
calculation to be described in detail. Additionally, the shortened or reduced 

20 filter updates may be implemented in an on-line fashion such that the 
reduced filter models can be determined without shutdown and without off- 
line pte-training. Apparatus are described herein for reducing the filter 
models in an on-line fashion for both MIMO systems and SISO systems, as 
well as for feed-forward and feed-back type systems. Further, the on-hne 

25 reduction method and apparatus may be used for feed-forward systems that 
utilize a synchronization signal and a waveform generator. 
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The above-mentioned and further novel aspects and advantages of 
the present invention will become apparent from the accompanying 
descriptions of the preferred embodiments and attached drawings. 

Brief Description of the Drawings 

5 The accompanying drawings which form a part of the specification, 

illustrate several embodiments of the present invention. The drawings and 
description together, serve to fully explain the invention. In the drawings: 

Fig- 1 is a block diagram illustration of a filtered-x active attenuation 
system known in the prior eirt; 

10 Fig. 2 is a schematic illustration of an active attenuation system with on- 
line feedback, error path, and output transducer modeling known in the 
prior art; 

Fig. 3 is a schematic illustration of an active attenuation apparatus and 
system of the present invention including an on-line reduction of the long 
15 filter model in a feed-forward system; 

Fig. 4A is a schematic illustration of an active attenuation system 
including on-line reduction of the filter model in a feed-back type control 
system (lacks a reference sensor); 

Fig. 4B is a schematic illustration of an active attenuation system including 
20 on-line reduction which includes a reduced filter model in a feed-forward 
control system which utilizes a waveform generator and synchronization 
signal; and 
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10 



15 



20 



Fig. 5 is a schematic illustration of an MIMO active attenuation system 
including on-line reduction for determining the reduced filter models in a 
feed-forward control system. 



In the following, it should be understood that the present invention 
may be used in systems for controlling acoustic vibrations in air as well as 
mechanical vibrations in gases, liquids, or solids. In summary, the 
present invention provides a method and apparatus for transforming a long 
filter model into a shortened or reduced filter model, thus reducing the 
computational burden required to update the control filters. The present 
invention is ideally suited, and is applicable to, tonal active control systems 
where one to four tones may be present and controlled. The term "tone" 
refers to the frequency of the mechanical vibration or noise. 

Fig. 3 illustrates one embodiment for carrying out the filter model 
reduction method. In this aspect of the present invention the reduction 
method is accomplished in an on-line fashion in reduction block 26. The 
active system 20 is used to cancel either a vibration or a sound and includes 
an update model 22 wherein the A Filter 32 and B-Filter 34 are adaptively 
updated to provide an output signal to an output device 25 in order to cancel 
a disturbance (the disturbing noise or vibration). Optionally, the system 20 
may also include on-line identification block 21 including on-line means for 
extracting an S'E' model 23. The S'E' model 23, i.e., the long filter model, is 
an estimate of the frequency response of the plant across the frequency 
range of interest. Basically, this is the long filter model Bk which will be 
described in the method to follow. This system 20 uses random low-level 
uncorrelated noise from random noise source 29 to develop the S'E' model 
23 on-line. The S'E' model 23 is a model of the error path 27 and the output 



Detailed Description of the Invention 
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device 25. However, it should be noted that the S'E' model 23 could be 



identification block 21. 

Once obtained, the S'E' model 23 is then copied and used in the 

5 novel on-line reduction process of the present invention taking place at 
reduction block 26. Within the reduction block 26, the input signal 36 is 
input into the reduced S'E' model 28 and to the S'E' copy 24. The input 
signal 36 in this embodiment is the signal that goes to both the A Filter 32 
and also to the reduction block 26. This convention for the input signal is 

10 the same for all of the embodiments to follow. The output from the S'E' 
model copy 24 is transferred through hne 38 and is summed (according to 
the sign convention shown) with the output from the reduced S'E' model & 
from line 40 at summing junction 42. The output from the summing 
junction 42 is input into convergence method 44 via line 46. Ideally, once 

15 the signal at 46 is zero or near zero, the reduced S'E' model 28 may be 
copied to reduced S'E' Copy blocks 30 and 30^ to be used in the update of the 
A Filter 32 and B Filter 34 within the update model 22. The error 
information 41 and the output from the reduced S'E' model copies 30 and 30' 
are used for update of the A Filter 32 and B Filter 34 models. It should be 

20 noted that the update calculation method which calculates the coefficients 
or weights for the A filter 32 and B filter 34 at adaptation junctions 31 and S 
may be suspended to allow the reduction method of the reduction block 26 to 
converge and arrive at the reduced S'E' model 28. The reduced S'E' model 
28 contains the exact same frequency response as the S'E' model 23, but 

25 only at the L number of discreet frequencies of interest, i.e., the tones to be 
controlled. 



optionally developed off-line by pre-training allowing elimination of 
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In the Fig. 3 embodiment, the apparatus obtains the input signal 36 



from a reference sensor (not shown) for sensing a source of input 
disturbance and providing an input signal indicative of the input 
disturbance. The input disturbance causes at least one of a vibration and 
5 noise to occur within a control volume (not shown). An output device 25 is 
used to provide a canceling output for locally canceling at least one of said 
vibration and noise at the defined location of the control volume. An error 
sensor (not shown) for sensing a residual combination of at least one of 
noise and vibration due to the canceling output and the input disturbance 

10 produces an error signal 41. An update model 22 is provided for 
transforming said input signal 36 into a canceling output which is 
transmitted to the output device 25. The control filter(s) 32 and 34 include a 
plurality of control filter weights to be updated. An x-filter S'E' model 23 is 
used for providing a frequency response estimate of the dynamic system or 

15 plant which is used in the reduction method. Finally, means are provided 
in reduction block 26 for shortening the x-filter S'E' model 23 such that the 
shortened or reduced filtered-x model 28 has a smaller nimnber of 
coefficients that match the frequency response of a plant at a finite number 
L of distinct disturbance frequencies. The smaller number of coefficients 

20 are used in an update-calculation method at junctions 31 and 33 which may 
incorporate a method such as filtered-x LMS, for adaptively updating the 
control filter weights of the control filter(s) 32 and 34. As a result of the 
fewer calculations needed, the updates to the A Filter 32 and B Filter 34 may 
be made more often or smaller CPU's may be utilized. 

25 The computational reduction is accomplished by a reduction 

method taking place in reduction block 26 whereby first an estimate of the 
plant or frequency response is obtained, i.e., the S'E' model 23. This 



estimate is referred to in this method as the long x-filter model Bk- This 
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can be done by exposing the system to an impulse function, random noise, 
or sine sweep and recording the response in an off-line or on-line fashion. 
The off-line procedure is known as pre-training aind generates a model of 
the frequency response of the plant over a wide frequency range. 

5 Further, the long x-filter model Bk may be initially obtained and 

also updated by on-line techniques, such as in on-line identification block 
21, where low level random uncorrelated noise from random noise source 
29 is input to the dynamic system and the frequency response is recorded. 
On-line updating is preferred for systems which the plant varies as a 
10 function of time, temperature, etc. Off-line pretraining may be appropriate 
where the system response characteristics are unchanging. For prior art 
filtered-x LMS systems, the output of the long x-filter model is then used 

directly in an update calculation method to calculate the updated feed- 
forward control filter weights. 

15 In the present invention, rather than using the output information 

rit from the long x-filter Bk directly, where the long filter model Bk has a 
length K, the long x-filter model Bk is transformed by a reduction method 
or technique into a reduced or shortened x-filter model Dm which has a 
length of M. In particular, in tonal systems where the method is 

20 preferable, there are a discreet number of dominant tones, where M<K and 
typically, M«K. 

The method of transformation will be illustrated with respect to a 
single input single output (SISO) system, but has equal applicability to 
MIMO systems where a distinct number of dominant tones are to be 
25 controlled. It should be noted that the advantage of the present invention 
method and apparatus will diminish as the number of frequencies of 
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interest increases. The output vector Yj consists of J consecutive samples 
from the long x-filter model Bk is described as 



yk yk-\ •■• yt-(y-i) 



where: 

5 yk~ individual sample output 

and Yj is obtained through on-line or off-line methods. Yj could 
then be used to calculate the shortened or reduced filter model 28 according 
to the following 

10 However, in practice, one would extract Bk first, which is an 

estimate of the plant. Then, an equivalent shortened or reduced x-filter 
model Dm would be constructed from the long x-filter model Bk- First, the 
matrix [HjjcJ would be generated by extracting a collection of inputs at 36. 
[Hjjj] represents a Hankel matrix which is a collection of row vectors, each 

15 row being a buffer of input data. In a feed forward system, the data may be 
collected from the reference sensor (not shown). In a feed-back type system, 
the input includes the error information. It should be understood that any 
appropriate input may be used to construct the Hankel Matrix. K 
represents the length of the long filter model Bk- Next, [Hj,m] would be 

20 extracted, as a subset of [HjjcJ- This is accomplished by taking the M left- 
most columns of the [Hj^ matrix. Next, one would complete or perform 
an inverse of [Hjjvi]- The inverse is denoted [Hj^+]. Methods of developing 
inverses are known, such as taught in Matrix Computations by Golub and 
Van Loan, Johns Hopkins University Press, 1983. 
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As mentioned above, the solution may be arrived at by calculating 
Yj as follows: 

where: 

- 

The solution to the equation above is the vector of coefficients Zk of 
the short x-filter Dm which can model the output response only at the L 
distinct, unique, and dominant tones present in the reference signal Xk- 

Two forms of the generaUzed solution which require a pseudo 
10 inverse operation are as follows: 
\st Form D^^ =[Hjj^ Hj jc]b^ 

2nd Form = [^j.m]^j 

These two forms are known as pseudo inverse methods where the 
generalized inverse (superscript +) is used to describe all solutions 
including cases where the actual inverse does not exist. When an actual 

15 solution does exist, for example, when [Hj^ is non-singular, the 
shortened x-filter model Dm and the Br long x-filter model will produce 
exactly the same output for the same input, but only at the L firequencies or 
dominant tones present in the reference signal sk- L is a user selected 
input that depends on the number of dominant tones present in the 

20 disturbance. 

For the first form, once [HjjvT^] and [Hjjd are computed, their 
product is multiplied by Br- Updates to [Hj^vi*] and [Hjjd are needed only 
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when there is a change to the system. The short model Dm contains far 
fewer than K terms present in the long x-filter model Bk for tonal systems. 
The first form is more desirable. The final step is to calculate the control 
filter weights for updating the A and B Filters 32 and 34, preferably 
5 according to a gradient descent update method based upon t>M and the 
error signal eu- In general, the ideal relationship for a tonal system is 
found to be where: 

2L<M<K 

The second form may also be used but is more computationally intensive. 

10 In particular, the present invention is a reduction method for 

reducing the computational burden of an active control system which uses 
any gradient descent control algorithm and which uses an estimate of the 
plant such as S'E' model 23. The reduction method has particular utility in 
systems using a filtered-x LMS algorithm. The method involves the steps of 

15 extracting an input signal xu at 36 , such as in a feed-forward system, 
indicative of characteristics of a source of disturbance. Next, an error 
signal 41 having magnitude eu indicative of a residual disturbance at a 
selected location is extracted. 

The reference signal is fed into a feed-forward control filter 32 for 
20 generating a control signal uk- The output signal uk is fed to an output 
device 25 to cancel the disturbance, that is, to minimize the magnitude ek of 
error signal 41. It is notable that the B Filter 34 may not be needed if the 
feed-back path 35 is not significant. The feed-forward control filter S 
within the update model 22 has a plurality of control filter weights which 
25 are adaptively updated according to an update calculation method taking 
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place at 31. The update-calculation method may be filtered-x LMS or the 
like. 

The transformation step taddng place at reduction block 26 involves 
transforming the long x-filter S'E' model 23 having a length of K into a 

5 shortened or reduced x-filter model 28 having a length of M, where M is 
less than K. Typically, the long filter model Br (the S'E' model 23) will 
model the plant over a wide frequency range, for example, from 100 Hz to 
400 Hz and may include 60-100 weights. The shortened x-filter model may 
only include 2-6 weights if the input only contains one or two tones to be 

10 canceled. Inputting the input signal 36 into said shortened x-filter model 28 
generates output 40. 

This output 40 is then used to converge the reduced filter model 2B 
by comparison with the output of the S'E' model 23 at the frequencies of 
interest. Then a copy 30 of the reduced filter model 28 is used to calculate 
15 the control filter weights of the feed-forward control filter 32 according to 
the update-calculation method at 31. The update-calculation method at 31 
uses both the shortened x-filter copy 30 output together with the error signal 
41. Finally, an output device 25 is driven responsive to the control signal uk 
in order to reduce at least one of vibration and sound at a selected location. 

20 The input signal at 36 in this embodiment is generated by a sensor 

such as at least one of a microphone, accelerometer, force transducer, and 
a tachometer, or the like. The error signal 41 in this embodiment may be 
generated by a sensor which can be at least one of a microphone, 
accelerometer, force transducer, or the like. Finally, the output signal is 

25 provided to an output device 25 that may be at least one of a loudspeaker, 
inertial actuator, active mounting, and electro-mechanical actuator, or the 
like. This reduction method has particular utility in controlling noise and 
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vibrations from sources having two dominant tones such as an aircraft 
turbofan engine. The first tone emanates from the low pressure turbine 
and the second from the high pressure turbine. The method is also useful 
for controlling 4 dominant tones in a turboprop aircraft. 

5 Fig. 4A illustrates an active system 2(y of the feed-back type. It is 

generally referred to as a feed-back type because there is no reference 
sensor present or needed in this embodiment. The system and apparatus 
utilizes the on-line reduction method of the present invention in the 
environment of a feed-back type system to obtain a reduced S'E' model S 

10 which is then used in the adaptation process at junctions 31 and 33. The 
on-line reduction is implemented in conjunction with the feed-back type 
system taught in US Pat. No 5,206,911 to Eriksson, the disclosure of which 
is hereby incorporated by reference herein. As in the feed-forward system 
of Fig. 3, the system 20* may include on-line identification block 21 for 

15 extracting an S'E' model 23 indicative of the error path 27 and the output 
device 25. Alternatively, the S'E* model 23 may be obtained with off-line pre- 
training and block 21 may be eliminated. In the Fig. 4A embodiment, the 
system 20' is identical to the system in Fig. 3 except that the Fig. 4A does not 
include a reference input, i.e., it is a feed-back type system. The input 

20 signal 39 for this embodiment is input into the A Filter 32 of the update 
model 22 and also to the reduction block 26. Note, the input signal 39 is 
derived from error signal 41' and the output from the reduced S'E* copy 30'. 

The on-line identification block 21, the update model 22, and the 
reduction block 26 are the same as in Fig. 3. Upon determining the S'E' 
25 model 23 and copying it to the S'E' Copy 24, convergence at junction 44 will 
result in a reduced S'E' model 28 which matches the frequency response of 
the S'E' model 23 at the distinct frequencies of interest. Therefore, the 
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reduced S'E' model 28 includes less information than the S'E' model 23. 
This is particularly true if the plant is dynamically dense or complex, for 
example, where the plant includes many resonances. The reduced S'E' 
model 28 is then copied to reduced copy models 30 and SC which are used 
5 for adaptively updating the coefficients of the A filter 32 and the B filter 34. 
The input 39 to the reduction block 26 is the combination of the error signal 
41' and the output of the reduced S'E' Copy 30'. The input signal 39 is 
directed to both the reduction block 26 and the update model 22. 



10 system and apparatus utilizing the on-line reduction technique of the 
present invention. In this particular embodiment, the feed-forw£u-d system 
includes a waveform generator 48 for generating at least one sinusoidal 
input signal at 40. The waveform generator 48 receives a synchronization 
signal 49 which contains the frequency content of the disttirbance. This 

15 synchronization signal 49 may include multiple dominant and discreet 
tones. By using on-Hne reduction block 26, the filter model S'E' 23, which 
contains frequency response information across the entire range, may be 
reduced in complexity and size. Therefore, updating of the coefficients of 
the A filter 32 is less computationally intensive. The method is essentially 

20 the same as the Fig. 4A and Fig. 3 embodiments except that the input 
signal 40 is the output of the waveform generator 48, 

Fig. 5 illustrates the present invention reduction method for 
transforming the long filter model into a reduced filter model is also 
applicable in a MIMO system 120. The description of this complex multi- 
25 channel system will reference only portions of the diagram, as much of the 
diagram is symmetric and duplicative. Like reference numerals have been 
utilized to signify like components as compared to the Figs. 3, 4A, and 4B 



Fig. 4B illustrates an active system 20" including a feed-forward 
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embodiments. In the Fig. 5 MIMO embodiment, each of the S'E' models is 
an estimate for the actual frequency response indicated as SEn, SE12, SE21, 
and SE22- For example, the SEn nomenclature used is indicative of the 
frequency response between the first output device si and the first error 
5 sensor ei. Further, the numeral 1 or 2 following RA indicates the sensor ri 
for subscript 1 and r2 for subscript 2. The numeral 1 or 2 following RB 
indicates actuator si for subscript 1 and S2 for subscript 2. Output device si 
may take the form of an actuator, speaker or other cancellation device. 

In a MIMO system of Fig. 5, the reduction method includes 
10 obtaining at least one input signal such as ri as is present in line 136 
indicative of characteristics of at least one source of disturbance such as an 
input noise. Further, at least one error signal such as e'l indicative of at 
least one residual disturbance (such as output noise) at least one location is 
extracted. The active system 120 also includes at least one output device 
15 such as si. The at least one input signal at 136 is input into at least one 
control filter 132 within an update model 122. The control filters are 
denoted as An, A21, A12 and A22, andBii,B2i, B12 and B22. The at least one 
control filter 132 generates at least one control signal. The at least one 
control filter 132 also has a plurality of control filter weights which are 
20 adaptively updated by an update-calculation method such as at 131. Again, 
the preferred method is filtered-x LMS or the like. 

At least one long filter model 123 such as S'E'n may be determined 
in an on-line fashion by on-line identification block 121 which represents at 
least one fi-equency response of said system. This may be also be 
25 accompUshed by off-line identification or pre-training. The at least one long 
filter model 123 which has a length of K is then transformed into at least 
one shortened or reduced filter model 128, such as RAiS'E'u which has a 
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length of M, where M is less than K. In particular, it is preferable that 
2L<M<K This is true of MIMO and SISO systems. 

The transformation for the S'E'n long filter model as it relates to ri 
takes place in reduction block 126. The reduction method used is identical 
5 to that described in relation to the Fig. 3 embodiment. Upon convergence of 
the reduction model in block 126, the reduced model is copied to the update 
model 122 to be used in the update process for updating the An control filter 
132. The reduction method is also performed for the remaunder of the 
frequency response models S'E'12, S'E'2i, and S'E'22 as they relate to ri. 
10 This occurs in upper A reduction Block 150 within upper reduction block 
152. The reduction of the B Filter models relating to actuator si is required 

in MIMO systems. These reductions take place in the Upper B reduction 
block 151 within upper reduction block 152 using input signal 157. 
Likewise, the reductions relating to sensor take place in Lower A 

15 reduction block 154 and Lower B reduction block 155 in lower reduction 
block 156. MIMO systems, which use the reduction techniques of the 
present invention are particularly useful for tonal cancellation systems 
where the multiple tones are to be canceled. By way of example and not by 
limitation, the system may include four or less dominant tones to be 

20 canceled. Upon generating the at least one shortened filter model 128, at 
least one input signal ri is input into said at least one shortened filter model 

128 to generate at least one shortened filter output. 

Next, the control filter weights of said at least one control filter 132 
are calculated according to said update-calculation method based upon said 
25 at least one shortened filter output and said at least one error signal. 
Finally, at least one output device, such as speaker si is driven responsive 
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to said at least one control signal in order to reduce at least one of vibration 
and sound at said at least one location, such as \\4thin the duct shown. 

It should be understood that this reduction method has applicability 
to MIMO systems where mechanical vibration is attempted to be controlled 
5 at a defined location by use of active mounts, actuators or the like. 

While several embodiments of the present invention have been 
described in detail, various modifications, alterations, changes and 
adaptations to the aforementioned may be made without departing from the 
spirit and scope of the present invention defined in the appended claims. It 
10 is intended that all such modifications, alterations, changes and 
adaptations be considered part of the present invention. 



21 



o 



o 



wo 96/14634 PCTAJS95yU453 

Claims 

What is claimed is: 

1. A method of reducing the computational burden of an active 
system, comprising the steps of: 

5 (a) obtaining an input signal indicative of characteristics of 

a source of disturbance, said input signal having dominant tones present 
therein to be canceled; 

Cb) extracting an error signal indicative of a residual 
disturbance at a selected location; 

10 (c) inputting said input signal into an update model having 

a control filter for generating a control signal, said control filter having a 
plurality of control filter weights which are adaptively updated according to 
an update calculation method; 

(d) determining coefficients of a long filter model which 
15 represents a fi-equency response of said system; 

(e) transforming said long filter model into a shortened 

filter model; 

(f) inputting said input signal into saud shortened filter 
model to generate a shortened filter output; 

20 (g) adaptively calculating said control filter weights of said 

control filter according to said update calculation method based upon said 
shortened filter output and said error signal; and 



22 



wo 96/14634 



PCT/US95/13453 



10 



(h) driving an output device responsive to said control 
signal in order to reduce at least one of vibration and sound at said selected 
location. 

2. A method of claim 1 wherein updating of said long filter model 
5 is accomplished in on-line fashion in an on-Une identification block. 

3. A method of claim 1 wherein said input signal in a feed 
forward system is provided by a reference sensor which is at least one of a 
microphone, accelerometer, force transducer, and tachometer. 

4. A method of claim 1 wherein said error signal is provided by 
an error sensor which is at least one of a microphone, accelerometer, and 
force transducer. 

5. A method of claim 1 wherein said output device is at least one 
of a loudspeaker, inertial actuator, active mounting, and electro- 
mechanical actuator. 

6. A method of claim 1 wherein said method is used to control 
two to four dominant tones in a tuxbofan aircraft and to control four to eight 
dominant tones in a turboprop aircraft. 

7. A method of claim 1 wherein said short filter model is obtained 
on-line in a reduction block. 

8. A method of claim 1 wherein said shortened filter model is 
calculated using a pseudo matrix inversion method. 

9. A method of claim 1 wherein a gradient descent method is 
used for providing an update calculation which is a filtered-x LMS 
algorithm. 
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10. A feed-forward control apparatus for active control, 
comprising: 

(a) a reference sensor for sensing a source of input 
disturbamce and providing an input signal indicative of said input 

5 disturbance, said input disturbance causing at least one of a vibration and 
noise to occur within a control volume; 

(b) an output device for providing a canceling output for 
locally canceling said at least one of said vibration and noise; 

(c) an error sensor for sensing a residual combination of at 
10 least one of noise and vibration due to said canceling output and said input 

disturbance and providing a residual error signal; 

(d) a control filter model for transforming said input signal 
into a canceling output, said control filter including a plurality of control 
filter weights; 

15 (e) a long filter model for providing an frequency response 

estimate; and 

(f) means for shortening said long filter model into a short 
filter model such that said short filter model has a smaller number of 
coefficients that match the frequency response of a plant at a finite number 
20 of disturbance frequencies, said smaller number of coefficients used in an 
update calculation for adaptively updating said control filter weights of said 
control filter. 

11. In an active control system utilizing a gradient descent 
algorithm, a method of updating control filter weights of a control filter, 

25 comprising the steps of: 
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(a) extracting a long filter model Bk which is an estimate of 

a plant; 

(b) generating a [Hj^k] matrix; 

(c) extracting a [Hj^mI matrix from said [Hj^kJ matrix; 

5 (d) completing an inverse of said [Hj^mJ matrix denoted as 

(e) calculating a shortened filter model Dm according to the 

relationship: 

Dm = [Hj^M^] * [Hj^ * Bk; 
10 and 

(f) calculating said control filter weights according to a 
gradient descent update method based upon an output fi'om said shortened 
filter model Dm and an error signal. 

12. A method of claim 11 wherein said calculation of said 
15 shortened filter model Dm is accomplished in on-line fashion in a reduction 

block. 

13. A method of claim 11 wherein said gradient descent update 
method for providing said update calculation is a filtered-x LMS algorithm. 

14. A method of claim 11 wherein said system is one of a feed- 
20 forward control system and a feed-back type control system. 

15. A method of claim 11 wherein said method is used in a feed- 
forward control system further including a waveform generator. 
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16. A method of reducing the computational burden of an active 
system, comprising the steps of: 

(a) extracting at least one input signal indicative of 
characteristics of at least one source of disturbance; 

5 (h) extracting at least one error signal indicative of at least 

one residual disturbance at least one location; 

(c) inputting said at least one input signal into at least one 
control filter for generating at least one control signal, said at least one 
control filter having a plurality of control filter weights which are 

10 adaptively updated according to an update calculation method; 

(d) determining coefficients of at least one long filter model 
which represents at least one frequency response of said system; 

(e) transforming said at least one long filter model into at 
least one shortened filter model; 

15 (f) inputting said at least one input signal into said at least 

one shortened filter model to generate at least one shortened filter output; 

(g) adaptively calculating said control filter weights of said 
at least one control filter according to sedd update calculation method based 
upon said at least one shortened filter output and said at least one error 

20 signal; and 

(h) driving at least one output device responsive to said at 
least one control signal in order to reduce at least one of vibration and 
sound at said at least one location. 
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17. A method of claim 16 wherein said method includes extracting 
a pliirality of input signals, extracting a plurality of error signals, and 
driving a pluradity of output devices in a MIMO system. 

18. A method of reducing the computational burden of an active 
5 system, comprising the steps of: 

(a) obtaining an input signal indicative of characteristics of 
a source of disturbance, said input signal having dominant tones present 
therein; 

(b) obtaining an error signal indicative of a residual 
10 disturbance at a selected location; 

(c) inputting said input signal into an update model having 
at least one control filter therein for generating a control signal, said at 
least one control filter having a plurality of control filter weights which are 
adaptively updated according to an update calculation method; 

15 (d) determining a long filter model which represents a 

frequency response of said system, said long filter model being obtained in 
one selected from an on-line and off-line fashion; 

(e) transforming said long filter model into a reduced filter 
model within a reduction block; 

20 (f) inputting said input signal into said reduced filter 

model to generate a reduced filter output within said update model; 

(g) adaptively calculating said control filter weights of said 
control filter within said update model according to said update calculation 
method based upon said shortened filter output and said error signal; and 
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(h) driving an output device responsive to said control 



signal in order to reduce at least one of vibration and sound. 

19. A method of claim 18 further including using said output 
signal as an input to a B reduction block for reducing said long filter model 

5 within said B reduction block to produce a reduced filter model to be used in 
updating a B filter within a MIMO system. 

20. A method of claim 18 further including using said input 
signal as an input to an A reduction block for reducing said long filter 
model within said A reduction block to produce a reduced filter model to be 

10 used in updating an A filter and using said output signal as an input to a B 
reduction block for reducing said long filter model within said B reduction 
block to produce a reduced filter model to be used in updating a B filter in a 
MIMO system. 
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